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INTRODUCTION  AND  REPORT  SUMMARY 


As  mankind  sets  its  sights  on  the  next  threshold  of  space 
exploration  (manned  expeditions  to  neighboring  planets)  ,  propulsion 
concepts  need  to  be  developed  to  meet  the  requirements  of  such 
missions.  Chemical  propulsion  systems  can  meet  the  criteria  needed 
for  nearby  planetary  trips,  but  nuclear  fusion  offers  better 
performance  in  specific  impulse  and  trip  duration  [Ref.  1] .  One 
fusion  propulsion  concept  which  the  Air  Force  is  studying  uses  the 
Dense  Plasma  Focus  (DPF) . 

The  DPF  propulsion  concept  was  first  studied  in-house  by  the 
Air  Force  at  the  Phillips  Laboratory  by  C.  Leakeas  [Ref.  1] .  His 
report  modeled  a  DPF  propulsion  system  in  its  entirety.  The  work 
herein  serves  to  streamline  and  update  the  previous  computer  model 
and  also  to  determine  the  effects  of  electromagnetic  (EM)  radiation 
losses  on  the  DPF  performance.  After  beginning  with  a  short 
discussion  of  the  cause  of  such  losses,  the  two  major  forms  of  EM 
radiation  losses  in  fusion  are  explained  in  detail,  those  being 
bremsstrahlung  and  synchrotron  radiation. 

A  primary  goal  of  this  report  is  to  investigate  the  effects  of 
the  synchrotron  radiation  losses  on  the  DPF  concept.  Special 
emphasis  is  placed  on  determining  the  synchrotron  radiation's 
emission  in  relation  to  the  plasma  pressure  ratio  parameter  (3.  The 
quantity  of  synchrotron  radiation  in  relation  to  bremsstrahlung  is 
mentioned  as  well.  Further  study  is  done  to  find  the  degree  to 
which  emitted  synchrotron  radiation  may  be  reabsorbed  by  the 
plasma . 

After  looking  at  the  effects  of  the  EM  losses,  the  remainder 
of  the  report  is  devoted  to  describing  the  new  DPF  computer  code. 
The  ideas  and  assumptions  used  in  the  original  program  are  given, 
along  with  the  alterations  and  reasons  for  their  presence  in  the 
new  code.  The  effects  on  DPF  rocket  performance  as  a  result  of  the 
new  code  are  mentioned  briefly.  Conclusions  of  the  current 
research  and  recommendations  for  future  research  in  both  the  areas 
of  DPF  EM  losses  and  DPF  computer  modeling  are  given  in  closing 
this  report. 
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GLOSSARY  OF  SYMBOLS 


a  =  Acceleration  [m-  s'2] 

ar  =  Radial  acceleration  [m-  s'2] 

B  =  Magnetic  flux  density  [T] 

Bi  =  Magnetic  flux  density  within  the  plasma  [T] 

B0  =  Outer,  or  surface,  magnetic  flux  density  [T] 

}  =  Plasma  pressure  ratio  parameter 

}a  =  Electron  pressure  ratio  parameter 
=  Ion  pressure  ratio  parameter 
c  =  Speed  of  light  in  a  vacuum  (2.99792458  [m- s'1] ) 
c'  =  Numerical  coefficient  [T~2-  s"1] 
cp,d  =  Specific  heat  of  deuterium  [ J-  kg'1*  “K'1] 
y  ’  =  Relativistic  mass  correction  factor 

Dp  =  Plasma  depth  [m] 

At  =  Length  of  time  in  one  pulse  [s] 

AT  =  Change  in  temperature  of  electrodes  during  one 

pulse  [K] 

e  =  One  electronic  charge  (1 . 6021892xl0'19  [C] ) 

E  =  Electric  field  strength  [V-m"1] 

Ee0  =  Electron  rest  mass  energy  (511.0  [keV] ) 

T|  =  Bremssthrahlung  classical  treatment  parameter  [Ref.  4] 
T\.  =  Electrical  resistivity  [C2-m] 

T\tot  =  Total  particle  collection  efficiency  in  plasma  pinch 
f (p)  =  Relativistic  Maxwellian  momentum  distribution  function 
F  =  Force  [N] 

g  =  Gaunt  relativistic  correction  factor  for 

Bremsstrahlung 

h  =  Planck's  constant  ( 6 . 626176xl0'34  [J- s]  ) 

Imax  =  Electrical  current  needed  to  create  required  magnetic 
flux  density  in  the  DPF  [A] 
k  =  Boltzmann's  constant  (1 . 380662xl0'23  [J-K-1]) 

K  =  Plasma  temperature  in  [keV] ,  (see  also  T  [°K] ) 

Ke  =  Electron  temperature  in  [keV] 

Ki  =  Initial  kinetic  energy  in  [keV] 

Kf  =  Final  kinetic  energy  in  [keV] 

K2()  =  Modified  Bessel  function  of  2nd  order 

K3()  =  Modified  Bessel  function  of  3rd  order 

X0  =  Characteristic  synchrotron  radiation  absorption 
length  [m] 

mav  =  Average  mass  of  a  particle  in  the  pinch  [kg] 

m3  =  Mass  flow  of  deuterium  [kg-  s'1] 

me  =  Rest  mass  of  electron  (9 . 1095346xl0'31  [kg]) 

^•0  =  Permeability  of  free  space  (4 n  x  10"7  [H-m'1]) 

nbeg  =  Initial  pinch  ion  density  [m'3] 

nD  =  Number  density  of  deuterons  [m'3] 

ne  =  Electron  number  density  [m-3] 

nt  =  Species  'i'  ion  number  density  [m'3] 

nHa  =  Number  density  of  3He  ions  [m'3] 

Ne  =  Number  of  electrons 

v  =  Photon  frequency  [s"1] 


2 


p 

Pr 

Pe 

Pz 

Pa 

p 

*-  brems 


*-  sync 
^tot 

p 

q 

q8 


Plnit 

T 

T 

e 

V 

V9 

Ve 

Vi 

Vlnic 

Vpln 

Z 

Zi 

Zi 

z2 

to 


Linear  momentum  [kg-m-s"1] 

Linear  momentum  in  radial  direction  [kg-m-s'1] 

Linear  momentum  in  tangential  direction  [kg-m-s'1] 
Linear  momentum  in  z  direction  [kg-m-s'1] 
Bremsstrahlung  power  loss  density  [W-m-3] 
Bremsstrahlung  power  [W] 

Cyclotron  power  [W] 

Ohmic  heating  power  [W] 

Radial  power  [W] 

Synchrotron  power  (=  dWe/dt)  [W] 

Total  radiated  power  (=  dW/dt)  [W] 

Particle  kinetic  pressure  [Pa] 

Total  electronic  charge  on  a  particle  [C] 

Electronic  charge  on  an  electron  (=  -e)  [C] 

Radius  of  gyration  [m] 

Plasma  radius  [m] 

Effective  reflectivity  of  emitted  synchrotron 
radiation 

Initial  mass  density  of  fuel  particles  in  DPF  entrance 
[kg-  m'3] 

Plasma  temperature  in  [°K],  (see  also  K  [keV] ) 
Electron  temperature  in  [°K] 

Velocity  [m-  s'1] 

Tangential  velocity  [m-  s'1} 

Speed  of  electron  [m-  s'1] 

Speed  of  particle  1  [m-  s'1] 

Initial  volume  of  rundown  (interelectrode)  chamber  [m3] 
Volume  of  plasma  pinch  [m3] 

Average  atomic  number  of  the  fusion  species 
Number  of  charges  on  ion  species  ' i' 

Number  of  charges  on  particle  1 
Number  of  charges  on  particle  2 
Angular  frequency  [s'1] 

Cyclotron  angular  frequency  [s'1] 

Electron  cyclotron  angular  frequency  [s'1] 

Electron  plasma  angular  frequency  [s'1] 
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ELECTROMAGNETIC  RADIATION  LOSSES  IN  THE  DPF 


ROLE  OF  ELECTRIC  AND  MAGNETIC  FIELD  LINES 

In  this  section,  the  fields  and  interactions  in  and  around  the 
DPF  plasma  pinch  region  are  analyzed.  Following  the  lead  of 
earlier  work,  the  pinch  will  be  approximated  as  having  a 
cylindrical  shape.  An  illustration  of  the  DPF  and  pinch  may  be 
found  in  Figure  1.  The  view  is  a  cross  section,  showing  the 
annular,  cylindrical  cathode  with  the  rod-like  anode  in  the  center. 
The  enclosed  portion  of  the  device  extending  from  the  end  of  the 
cathode  to  the  nozzle  entrance  is  the  mixing  chamber.  In  it,  the 
remaining  fusion  reactants  and  products  mix  with  the  propel lant, 
heating  it.  This  results  in  a  high  temperature  gas  mixture,  which 
is  expanded  through  the  nozzle  to  produce  thrust.  A  dash-lined 
ellipse  surrounds  the  area  known  as  the  pinch.  When  necessary,  the 
pinch  dimensions  used  in  this  report  are  radius  =  0.0015  [m]  and 
length  =  0.0254  [m]  .  These  are  the  values  assumed  by  Leakeas  [Ref. 
1]  • 


Charged  particles  exhibit  changes  in  their  motions  and  emit 
radiation  due  to  additional  kinetic  energy  imparted  through 
accelerations.  These  accelerations  arise  from  interactions  with 
electromagnetic  fields.  Because  plasmas  are  composed  of  ions  and 
electrons,  interactions  involving  them  are  the  key  concerns.  The 
most  prominent  examples  are  the  interaction  of  electrons  with 
electric  fields  (bremsstrahlung)  and  the  interaction  of  electrons 
with  magnetic  fields  (synchrotron  radiation) .  Both  of  these 
radiation  types  will  be  covered,  with  a  special  emphasis  given  to 
synchrotron  radiation,  defined  in  the  SYNCHROTRON  RADIATION 
section.  Its  dependence  on  the  parameter  (5  for  emission  is 
detailed  in  the  EFFECT  OF  PLASMA  BETA  ON  SYNCHROTRON  RADIATION 
section.  Its  reabsorption  following  emission  is  discussed  in  the 
ABSORPTION  OF  SYNCHROTRON  RADIATION  section. 


BREMSSTRAHLUNG 

Bremsstrahlung,  German  for  "braking  radiation"  [Ref.  2],  is 
the  radiation  emitted  by  a  moving  charged  particle  accelerated  by 
collisions  with  other  particles.  Because  a  particle's  acceleration 
is  inversely  proportional  to  its  mass  for  a  given  force,  the 
bremsstrahlung  emitted  by  ions  is  much  less  than  that  of  electrons. 
The  electrons,  however,  are  deflected  (accelerated)  much  more  by 
the  ions  present  than  by  other  electrons.  Thus,  in  the  field  of 
fusion,  bremsstrahlung  is  usually  considered  to  emanate  from  the 
interaction  of  moving  electrons  with  essentially  stationary  ions. 

Moving  electrons  constantly  come  in  contact  with  the  electric 
fields  arising  from  other  charged  particles  (either  ions  or  other 
electrons)  .  Bremsstrahlung  can  be  emitted  through  a  process  known 
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Figure  1 

Dense  Plasma  Focus  Overview 


as  a  "free-free"  transition  of  an  electron.  In  this  process,  a 
free  electron  is  deflected  by  Coulomb  forces  due  to  other  charged 
particles  in  its  vicinity  [Ref.  3]  .  The  energy  imparted  to  the 
electron  during  the  deflection  is  released  in  the  form  of  a  photon, 
whose  energy  is  governed  by: 


K±  -  Kg  +  hv 


(1) 


with  Kx  equal  to  the  initial  kinetic  energy  of  the  electron,  Kf 
equal  to  its  final  kinetic  energy,  and  hv  equal  to  the  energy 
released  in  the  photon. 

When  the  energy  released  in  the  photon  approaches  the  kinetic 
energy  of  the  emitting  electron,  a  quantum  treatment  is  necessary 
to  gain  a  full  understanding  of  the  process  [Ref.  3]  .  If  being 
relativistic  is  considered  as  having  a  speed  one-tenth  that  of 
light  or  greater,  the  energy  at  which  electrons  attain  such  status 
may  be  found.  An  electron's  rest  mass  (E,0)  of  511.0  [keV]  is  equal 
to  mec2 .  Its  kinetic  energy,  accounting  for  a  relativistic  increase 
in  mass,  is  (y-l)mec2.  The  equation  for  the  relativistic  correction 
factor  y  is: 


Y 


_ 1 

^1  -  (vl/c2) 


Ea0 


(2) 


For  v„  =  0.1c,  y  =  1.005038,  with  K„  =  2.574  [keV]  by  Eq.  (2). 
Because  this  threshold  is  below  the  temperature  currently  sought 
for  operating  the  DPF,  one  must  consider  a  relativistic  treatment 
of  the  bremsstrahlung  energy  loss. 

Further,  Jackson  [Ref.  4]  states  that  a  classical  treatment 
of  bremsstrahlung  is  only  valid  if: 


T) 


(10*7) 


Z1Z2e2c2 


>  1 


(3) 


where  Z1  is  the  number  of  charges  on  the  first  particle  (an 
electron,  so  Zx  =  1),  Z2  is  the  number  of  charges  on  the  second 
particle  (a  3He  ion,  so  Z2  =  2),  e  =  1.6022xl0'19  [C],  h  is  Planck's 
constant  ( 6 . 6262xl0-34  [J-s] )  ,  and  vx  is  the  velocity  of  the  first 
particle.  For  an  electron  of  kinetic  energy  2.574  [keV]  (ve  = 
0.1c)  as  mentioned  above,  the  corresponding  speed  of  the  electron 
is  2.9979xl07  [m/s].  Thus,  T]  =  0.0232,  which  shows  that  a 

classical  treatment  is  invalid,  confirming  the  previous 
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relativistic  assumption. 

One  will  find,  however,  that  a  classical  approach  gives  the 
correct  functional  dependence  for  most  of  the  physical  parameters 
[Ref.  3]  .  A  quantum  correction  is  often  handled  by  a  correctional 
coefficient  known  as  the  Gaunt  factor,  denoted  by  '  g'  .  This  factor 
depends  on  the  type  of  particles  interacting.  In  a  fusion  plasma, 
the  main  interaction  is  that  of  electrons  with  the  Coulomb  fields 
of  charged  nuclei.  For  it,  a  value  of  g  =  1.2  is  commonly  used. 
The  bremsstrahlung  power  density  [W/m3] ,  including  Gaunt  factor,  is 
given  to  be  [Ref.  3]  : 


PB  -  5 . 35xlO"37^/Tg  ng  EnjZ] 


(4) 


where  Ka  is  the  electron  temperature  in  [keV] ,  ne  is  the  number 
density  of  electrons  in  the  plasma  [m~3] ,  nt  is  the  number  density 
of  ion  species  '  i'  in  the  plasma  [m“3],  and  ZL  is  the  atomic  number 
of  ion  species  'i'. 


SYNCHROTRON  RADIATION 

Synchrotron  radiation  originates  from  the  acceleration  of 
charged  particles  as  a  result  of  the  non-alignment  of  their 
velocities  with  their  local  magnetic  fields.  A  derivation  of  the 
expression  for  synchrotron  radiation  begins  with  the  equation  for 
the  force  acting  on  the  charge  in  question  due  to  electric  and 
magnetic  fields: 


F  -  q[E  +  vxB) 


(5) 


Considering  no  Couloumb  field  E  to  be  present  and  the  force  F  to  be 
given  by  mass  times  acceleration,  with  centripetal  acceleration 
equal  to  the  tangential  velocity  (v0)  squared  over  the  radius  of 
gyration  -.out  the  B  field  line,  the  resulting  form  of  Eq.  (5)  is: 


F 


qv9B 


(6) 
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Vector  notation  has  been  discarded,  as  va  is  orthogonal  to  B.  As 
stated  by  Tamor  [Ref.  5]  and  as  will  be  seen  later,  the  synchrotron 
loss  rate  varies  inversely  as  the  cube  of  the  particle  mass.  Thus, 
the  contribution  of  the  synchrotron  radiation  arising  from  ions  may 
be  neglected  in  comparison  with  the  contribution  from  electrons. 
The  frequency  of  emission,  the  cyclotron  frequency,  is: 


0) 


QB 

ym 


(7) 


and  for  an  electron,  q  =  e,  m  =  me,  and  (0C  =  G)ce. 

For  a  charge  q  (here,  an  electron,  qe)  moving  with  velocity  v 
and  acceleration  a,  the  total  radiated  power  is  given  by: 

p  m  dw  _  q2  a2  -  ( vxa)  /  c2 
tot  dt  6ne0c  (1  _  j/^3  '  (8) 

c2 


where  |a|  =  |ar|  and  |a|  =  |o)cav0|  =  |  (qaB/Yma)  v0|  [Ref.  3].  The 

diagram  in  Figure  2  shows  the  locations  of  the  vectors  in  the 
cylindrical  (r-0-z)  coordinate  system  being  employed,  assuming  a 
nearly  straight  B  field  line.  The  electrons  spiral  about  the  B 
field  lines,  tracing  out  elliptical  paths.  Eq.  (8)  becomes: 


P 


tot 


Qp  B2  Vg2 

6  7i  e0  cy2m2 


(9) 


With  radial  momentum  pr  =  7rnev0,  the  expression  for  total  radiation 
from  an  electron  without  consideration  of  whether  it  escapes 
follows  from  Eq.  (9)  to  be: 


P 


sync 


dwg 

dt 


do  B2p  2 
6neocy*m0* 


(10) 
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Now,  (pr2/2(ome)  =  NekTe  and  for  an  isotropic  distribution,  pr2 
=  (2/3)<p2>.  If  (dWe/dt)f(p)  is  now  computed  for  Eq.  (10)  with  f  (p) 
being  the  ralativistic  Maxwellian  distribution  function,  Rose  and 
Clark  [Ref.  6]  have  shown  that: 


Ho  e*B2  NakTa  K3(mac2/kTe) 
3nma2c  nta  K^lm^/kTj 


(11) 


where  K2  ( . . . )  and  K3(...)  are  Modified  Bessel  functions  of  second 
and  third  order,  respectively.  The  inverse  dependency  of  the 
synchrotron  radiation  energy  loss  on  the  cube  of  the  mass  is  now 
clearly  evident,  showing  the  reason  for  the  neglect  of  the  losses 
due  to  ion  interactions. 

Using  the  electron  cyclotron  radial  frequency  (Ci)ce)  defined 
earlier,  the  equation  for  synchrotron  power  emission  becomes: 


sync 


Ho  e2oca  (  Na  kTa )  (i  + 


Kn 


371 C 


JR. 


204.5 


+  .  .  . ) 


(12) 


of  which  the  cyclotron  radiation  power  (Pcyc)  portion  comprises  only 
the  first  term,  that  is: 


P 


eye 


Hq  Na  kT0 

3tc  c  ma 


(13) 


The  units  for  Ke  are  [keV] ,  while  the  units  for  the  other  terms  may 
be  found  in  the  GLOSSARY  OF  SYMBOLS.  If  one  assumes  that 
relativistic  effects  become  important  when  the  value  of  Psync  exceeds 
the  value  of  Pcyc  by  10%,  then  the  Te  value  at  which  this  occurs  is 
found  by: 


P  ■  1  1  P  m  P  ( 

^  sync  •  ^  eye  ^  eye ' 


i  + 


204.5 


) 


(14) 


which  yields  Ke  «  20  keV.  The  temperature  regime  being  considered 
for  the  DPF  model  in  this  report  is  about  50  to  100  keV.  One  sees, 
then,  that  a  relativistic  treatment  of  this  radiation  type  is 
required. 

The  nomenclature  associated  with  synchrotron  radiation  has 
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become  ambiguous.  Technically,  cyclotron  and  synchrotron  radiation 
arise  from  the  same  mechanism,  but  their  values  can  be  quite 
different.  Stated  simply,  cyclotron  radiation  is  the  non- 
relativistic  form  of  the  expression  for  synchrotron  radiation 
emission;  that  is,  the  classical  treatment.  For  cyclotron 
radiation,  the  frequency  of  emission  is  equal  to  the  gyrational 
(cyclotron)  frequency  of  the  plasma.  Synchrotron  radiation,  on  the 
other  hand,  can  be  emitted  at  frequencies  many  times  that  of  the 
plasma  cyclotron  frequency.  Because  the  operating  temperature 
regime  of  the  DPF  model  used  in  this  report  is  above  the 
relativistic  threshold,  cyclotron  radiation  is  not  an  appropriate 
name  for  this  type  of  radiation.  Thus,  the  term  synchrotron 
radiation  will  be  used  in  the  remainder  of  this  report,  as  it 
correctly  applies  to  all  energy  levels  of  this  radiation  type,  much 
as  relativistic  mechanics  can  be  used  to  describe  classical 
mechanical  situations. 


EFFECT  OF  PLASMA.  BETA  ON  SYNCHROTRON  RADIATION  EMISSION 

To  understand  the  role  of  synchrotron  radiation  in  modeling 
the  DPF,  the  parameter  (3  shall  now  be  introduced.  The  P  value  of 
a  plasma  is  the  ratio  of  the  particle  kinetic  pressure  to  the 
confining  magnetic  pressure  at  the  boundary.  This  is  given  by: 


P 


P 

J9j/2|X0 


(15) 


where  is  the  permeability  of  free  space  (4it  x  10'7  [H/m] )  ,  P  is 
the  particle  kinetic  pressure  [N/m2] ,  and  B0  is  the  magnetic 
induction  [T]  at  the  plasma  surface.  When  (3  is  equal  to  1,  a 
pressure  balance  between  the  expanding  plasma  particle  pressure  and 
the  confining  magnetic  pressure  is  in  existence.  A  value  of  p  less 
than  1  indicates  that  the  confining  magnetic  pressure  must  be  made 
stronger  for  a  given  particle  pressure  in  order  to  prevent  pinch 
expansion.  To  re-iterate,  it  is  assumed  that  a  uniform 
constriction  is  present  along  the  length  (0.0254  [m] )  of  the  pinch, 
giving  it  a  cylindrical  shape. 

A  P  value  of  1  requires  the  minimum  magnetic  induction 
necessary  to  stablize  the  plasma  for  a  given  value  of  P.  (3  =  1  is 
a  common  assumption  when  dealing  with  the  DPF,  as  modeling  has 
shown  such  a  p  value  to  be  stable,  even  in  the  presence  of  slight 
oscillations  [Ref.  7]  .  The  beneficial  consequence  of  such  a  choice 
for  P  is  the  elimination  of  losses  due  to  synchrotron  radiation. 
Hulme  [Ref.  8]  states  that  for  electrons  having  a  Maxwellian 
distribution,  the  power  density  due  to  synchrotron  radiation  (Psync) 
is : 
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(16) 


sync 


-  rJ 


C'  *  Bx  •  n 


where  c'  is  a  numerical  coefficient  with  units  [T'2*  s*1] ,  Bi  is  the 
magnetic  flux  density  within  the  plasma  (assumed  constant) ,  n,  is 
the  number  of  electrons  per  unit  volume,  and  T,  is  the  temperature 
of  the  electrons. 

Now  assume  that  both  the  ions  and  electrons  are  at  a 
temperature  T.  For  a  50-50  D-3He  mixture,  the  number  of  electrons 
is  equal  to  1.5^,  where  nL  is  the  number  of  Deuterium  (nD)  and 
Helium-3  (nHe)  ions  in  the  plasma:  n£  =  nD  +  nHe.  The  neutral 

plasma  pressure  equation  may  be  written  as  [Ref.  2] : 

P  -  ^'{nm+n±)  -\kT  -  (fl.+n*)  kT  ,  (17) 


where  n„  is  the  electron  particle  density  [m*3] ,  nt  is  the  ion 
particle  density  [m*3] ,  and  (3/2) kT  is  the  average  charged  particle 
three-dimensional  energy  [J] .  If  we  now  equate  the  plasma  pressure 
to  the  pressure  exerted  on  the  plasma  by  the  containing  (outer) 
magnetic  field,  B0,  Eq.  (17)  becomes: 


n2 

(n^+n^  kT  -  — -2- 

*  1  2ji0 


(18) 


In  addition  a  term  may  be  added  to  the  plasma  pressure  to 
account  for  any  residual  magnetic  fields  within  the  plasma.  This 
term  is  based  on  the  inner  magnetic  flux  density  and,  when 
inserted,  Eq.  (18)  becomes: 

(n.*nt)  kT  -  B\"  B‘  .  (19) 


Letting  (3  =  1—  BL2/B02,  we  have  the  expression  [Ref.  8]: 

(nQ  +  ni)  kT  -  ,  (20) 

2n0 
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which  may  be  compared  to  our  original  definition  of  0  in  Eq.  (15) . 
Eq.  (18)  can  be  obtained  by  letting  (5  =  1,  the  condition  in  which 
no  residual  magnetic  fields  remain  within  the  plasma.  In  terms  of 
Eq.  (20)  becomes: 


{n^+n^kT  - 


2Ho 


(21) 


This  expression  may  be  solved  for  Bj2,  with  the  resulting  expression 
substituted  into  Eq.  (16)  above  to  obtain: 


Psync  -  c'  '  *  Ini+IlJlJ  T2 


(22) 


It  is  now  apparent  what  a  P  value  of  1  means.  The  right  hand 
side  of  Eq.  (22)  becomes  0,  which  means  the  synchrotron  losses 
escaping  the  plasma  fall  to  0.  This  is  an  important  result,  for 
the  escaping  energy  due  to  synchrotron  radiation  can  be  quite  large 
and,  thus,  quite  detrimental  to  a  fusion  reactor's  ability  to 
ignite  and  continue  operating.  The  effect  of  (5  on  power  loss  due 
to  synchrotron  radiation  is  further  confirmed  by  the  work  of  S. 
Tamor  [Ref.  5] . 

Tamor  presents  a  relation  for  the  ratio  of  synchrotron  power 
to  bremsstrahlung  power  [Ref.  5],  given  by: 


°-3^~r  4  '  (23) 

0 


P 

*3ync  — 
p 

*  hrema 


where : 


ni  -  0.01-73.  (y-^) 


z* 


‘-ft’’ 


(24) 
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K  is  the  temperature  of  the  plasma  [keV] ,  rp  is  the  radius  of  the 
plasma  [m] ,  R  is  the  effective  reflectivity  for  emitted  synchrotron 
radiation,  Z  is  the  average  atomic  number  of  the  fusion  fuel 
species,  and  px  and  P„  are  the  plasma  beta  values  for  the  ions  and 
electrons,  respectively.  One  sees  that  Eq.  (23)  is  dependent  on 
the  factor  (l~P)/p.  When  (3  =  1,  P,ync/pbrM>«  =  0*  Thus,  no 
synchrotron  radiation  is  emitted  for  0  =  1,  agreeing  with  the 
result  found  in  Hulme  [Ref.  7] . 


ABSORPTION  07  SYNCHROTRON  RADIATION 

It  is  assumed  in  this  report  that  the  value  of  R  is  0,  which 
is  a  worst  case  situtation.  In  order  to  see  the  effect  of  a  P 
value  less  than  1,  an  arbitrary  value  of  {3  =  0.9  will  be  assumed. 
A  mass  density  of  2.2xl0"4  [kg-irf3]  is  used  for  the  pinch  volume  of 
the  DPF  [Ref.  1].  Using  Eqs .  (23,24)  with  the  parameters:  R  =  0, 
L  =  1.5  [cm]  (pinch  radius),  (3  =  0.9,  Z  =  1.5  (equal  amounts  of  D 
and  3He)  ,  nD  =  n3Ha  =  1.179xl026  [m-3] ,  and  pe  =  1.50!,  Psy„c/Pbrems  values 
for  a  temperature  range  of  10  to  100  [keV]  are  calculated.  This  is 
the  approximate  range  of  validity  Tamor  gives  for  Eqs.  (23,24). 
These  values  are  presented  graphically  in  Figure  3 .  One  sees  that 
for  a  p  value  of  0.9,  the  synchrotron  power  escape  from  the  pinch 
remains  modest  (P3ync  <  0 . 1- Pbrams)  until  about  62  [keV]  .  As  the 
temperature  increases,  however,  reabsorption  becomes  increasingly 
important,  as  is  evident  in  Figure  3.  One  sees  that  the 
synchrotron  power  is  climbing  significantly  (as  temperature 13/4)  as 
the  bremsstrahlung  power  is  increasing  slowly  (as  temperature1'2)  . 

In  Figure  4  one  sees  further  the  effect  of  lowering  the  P 
value  in  the  device  to  0.9.  Assuming  the  synchrotron  radiation  is 
not  reabsorbed,  the  DPF  is  limited  to  an  operating  regime  of 
approximately  34  to  117  [keV] ,  the  range  in  which  the  charged 
particle  power  exceeds  the  total  radiative  loss  power.  Figure  4 
indicates  that  the  optimum  operating  temperature  would  be 
approximately  75  [keV] ,  at  which  point  the  charged  particle  power 
exceeds  the  radiative  loss  power  by  the  greatest  amount.  It  has 
been  assumed  that  the  equations  of  Tamor  [Ref.  5]  remain  valid  up 
to  a  plasma  temperature  of  130  [keV] ,  in  order  to  show  the  power 
curve  intersections.  Achieving  a  p  value  of  1  eliminates  the 
necessity  of  reabsorption  by  preventing  the  escape  of  synchrotron 
radiation.  It  also  allows  one  to  avoid  the  concern  for  reflective 
walls,  which  is  not  specifically  addressed  in  this  report.  In 
Figure  4  one  notices  that  from  about  85  keV  upward,  the  net  power 
gain  (the  difference  between  the  P_brems  and  P_charged_particles 
curves)  for  a  P  =  1  condition  remains  approximately  constant.  A  P 
=  1  condition  is  also  the  most  desirable  approach  in  achieving 
maximum  DPF  performance,  as  is  shown  in  the  next  section. 

Even  if  the  P  value  is  less  than  1,  the  plasma  absorption 
capability  is  usually  sufficient  to  achieve  the  99.9%  synchrotron 
radiation  reabsorption  acceptance  level  stated  by  Dawson  [Ref.  9] . 
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Figure  3 

P_brems  [MW],  P_sync/P_brems  vs.  K  [keV]  for  n_e  =  3.534E+26 


Power  [MW] 


K  [keV] 

Figure  4 

Power  [MW]  vs.  K  [keV]  for  n_e  =  3.534E+26  m 
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A  body  that  emits  radiation  is  also  an  absorber  of  radiation.  The 
mechanisms  of  emission  and  absorption  for  the  transfer  of  radiation 
are  associated  with  one  another  by  the  Kirchoff  relation  [Ref.  9] : 

end s si  vi  m  q2  (25) 
absorptivity  9 


This  is  valid  so  long  as  the  electrons  in  the  plasma  are  described 
by  a  Maxwellian  distribution.  Here,  this  has  been  assumed  to 
extend  to  the  relativistic  Maxwellian  distribution  assumption 
stated  earlier.  Further,  the  expression  for  the  characteristic 
absorption  length  is  given  as  [Ref.  5]: 


c  (ii 


to 


1 .7  xlO18 


Bo 


(26) 


By  substituting  into  Eq.  (20)  and  considering  n1=(2/3)ne  for  a  50-50 
D-3He  fuel  mixture: 


(1.7 xlO19)  \i0kT 
2BoX0 


(27) 


one  sees  that  P  is  inversely  proportional  to  X0.  Thus,  a  decrease 
in  the  plasma  beta  results  in  an  increase  in  the  characteristic 
reabsorption  length. 

For  a  typical  50-50  D-3He  pinch  with  P  =  1,  ne  =  3.534xl026  [m' 
3] ,  and  optimum  operating  temperature  K„  =  75  [keV] ,  a  B0  value  of 
4217  [T]  is  obtained,  which  yields  X0  =  2.138xl0~5  [m]  .  The  reader 
is  reminded  that  a  P  =  1  condition  means  Bi  =  0,  as  shown  earlier. 
A  typical  radius  assumed  for  the  plasma  pinch  is  1.5xl0-3  [m] ,  which 
is  70.2  times  greater  than  X,0,  showing  why  all  of  the  synchrotron 
radiation  is  considered  to  be  reabsorbed  within  the  confines  of  the 
plasma.  This  confirms  the  earlier  results  found  using  the 
equations  in  Hulme  [Ref.  7]  and  Tamor  [Ref.  5]  and  presents  the 
physical  basis  for  P  =  1  leading  to  Psync  =  0.  Assuming  the  size  of 
the  pinch  to  remain  relatively  stable  to  changes  in  plasma 
temperature,  magnetic  induction,  and  particle  density,  the  K„  value 
at  which  X0  equals  the  plasma  pinch  radius  may  be  found.  This 
occurs  when  B/n„  =  8.824xl0*22  [T- m3]  .  By  Eq.  (18),  with  n*  =  (2/3)ne 
for  a  50-50  D-3He  plasma,  the  corresponding  electron  temperature  is 
based  on  the  relation: 
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(28) 


Ke  -  1 . 315  B0  . 


For  the  B0  value  of  4217  [TJ ,  K„  =  5545  keV,  or  5.545  MeV.  This 
value  is  much  larger  than  the  operational  temperatures  being  sought 
after  for  use  in  the  DPF  [Ref.  10]  .  This  indicates  that  the  escape 
of  synchrotron  radiation  from  the  dense  plasma  focus  device 
operating  with  J3  =  1  is  not  a  major  concern  and  is  justifiably 
neglected. 
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DENSE  PLASMA  FOCUS  COMPUTER  MODEL 


THE  DENSE  PLASMA  FOCUS:  A  FUSION  PROPULSION  SCENARIO 

The  recent  work  of  Chris  Leakeas  [Ref.  1]  has  laid  the 
groundwork  for  the  DPF  computer  model  used  in  this  study.  His  code 
was  entitled,  "The  Dense  Plasma  Focus:  A  Fusion  Propulsion 

Scenario."  The  newer  model  is  similar  in  some  respects  to  the 
previous  model.  The  older  code  was  based  on  a  DPF  system  utilizing 
various  pulsed  mode  operations  and  a  D-3He  fuel  mixture. 
Currently,  the  new  code  studies  the  nearly  steady  state  operation 
of  the  DPF,  albeit  still  in  a  pulsed  mode.  As  in  the  previous 
program,  the  mission  focus  was  determination  of  the  DPF's 
feasibility  for  a  trip  to  Mars. 

Leakeas'  code  was  written  with  a  linear  program  structure.  A 
number  of  key  assumptions  about  the  propulsion  device's  operation 
and  technological  feasibility  were  made.  These  are  carried  over 
into  the  current  model.  They  are: 

1)  Plasma  pinch  temperature  scales  as  current  squared. 

2)  Pinch  dimensions  are  length  0.0245  [m]  and  radius  0.0015  [m] . 

3)  Ions  reach  thermal  equilibrium  inside  the  pinch,  allowing  use 
of  Maxwellian  reaction  rate  parameters. 

4)  High  temperature  materials  beyond  those  currently  available 
will  be  developed. 

5)  Electrodes  and  mixing  chamber  (nozzle  entrance  area)  walls  can 
be  sufficiently  cooled  to  avoid  damage. 

6)  Propellant  becomes  completely  dissociated  and  ionized  in  the 
mixing  chamber  at  5000  K. 

7)  Capacitor  technology  advancement  will  allow  factor  of  10 
increases  in  specific  energy  over  present-day  capacitors  (0.2 
kJ/kg)  and  discharging  rates  of  100  Hz. 

8)  Confinement  times  can  be  increased  by  a  factor  of  100  to  about 
1 0~4  seconds. 

9)  Magnetic  fields  applied  downstream  (nozzle  entrance  area)  from 
the  pinch  do  not  adversely  affect  the  pinch  formation  or 
confinement  time. 

Many  of  the  findings  in  Leakeas'  report  were  quite  favorable 
toward  the  DPF's  use  in  a  Mars  mission.  Calculations  were  made 
within  the  framework  of  the  system  analysis  that  involved  the 
absorption  and  reflection  of  bremsstrahlung  and  synchrotron 
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radiations.  I3p  values  in  the  neighborhood  of  104  seconds  and 
initial  thrust-to-weight  ratios  in  the  proximity  of  10-3  were 
calculated.  The  lifted  payload  was  10s  kg. 


A  DENSE  PLASMA  FOCUS  MODEL  (IMPROVED  DPF  CODE) 

The  culmination  of  writing  the  new  DPF  code  was  the  FORTRAN 
program  "A  Dense  Plasma  Focus  Model."  Some  programming  changes 
were  made  in  creating  the  new  code.  For  instance,  many  of  the 
variable  names  are  now  mnemonic.  One  example  is  that  all  variables 
for  power  (energy  per  unit  time)  begin  with  the  letter  'P'  or  the 
letters  '  POW' .  Also,  all  variables  regarding  time  begin  with 
'TIM',  while  those  for  temperature  begin  with  'T'  alone.  The 
variable  descriptions  at  the  beginning  of  the  code  listing  in 
Appendix  A  should  clarify  any  further  questions  on  this  matter,  as 
every  variable  used  in  the  program  is  explained,  accompanied  by  its 
units . 


Program  execution  is  governed  by  a  set  of  SUBROUTINE'S  in  the 
main  program.  One  key  note  is  that  the  program  will  stop  running 
if,  after  initial  calculations,  it  finds  that  the  power  being 
generated  by  the  turbine/generator  system  is  not  sufficient  to 
recharge  the  capacitor  banks  fully.  This  is  a  critical 
requirement,  for  the  electrical  current  needed  by  the  DPF  for  each 
pulse  must  be  supplied  in  order  for  operation  of  the  DPF  to 
continue  at  the  wanted  keV  level  of  the  pinch. 

A  major  concern  which  arose  during  the  evolution  of  the 
current  code  was  the  fact  that  all  cases  studied  showed  the  need 
for  an  additional  power  source  to  recharge  the  capacitor  banks 
sufficiently.  The  reason  for  this  finding  is  the  result  of  the 
newer  treatment  given  many  aspects  of  the  DPF.  Some  differences 
were  minor,  but  others  were  quite  detrimental  to  the  performance  of 
the  DPF. 

One  area  of  change  was  the  calculation  of  the  mass  of  the 
system  used  to  store  and  pump  the  hydrogen  propellant.  The 
original  assumption  was  that  the  propellant  pumping  and  storage 
system  should  have  a  mass  equal  to  15%  of  the  propellant  mass. 
This  was  changed  to  10%,  the  figure  used  in  determination  of  the 
fusion  fuel  pumping  and  storage  system  mass,  based  on  the  mass  of 
the  fuel  (D  and  3He) .  This  would  also  decrease  the  mass  of  the 
DPF-powered  vehicle,  increasing  its  performance. 

Discussions  with  Dr.  C.  K.  Choi  of  Purdue  [Ref.  11]  led  to 
including  the  D-D  side  reactions  in  the  power  density  and 
electromagnetic  radiation  loss  calculations.  He  felt  that  for  the 
short  duration  (10-4  seconds)  of  the  pinch,  the  D-D  reactions  would 
be  significant,  but  secondary  reactions  such  as  D-T  would  not  be. 
Higher  power  densities  can  be  achieved  with  increased  amounts  of 
3He  in  the  fuel  mix,  but  at  a  cost  of  increased  fuel  mass  and 
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increased  bremsstrahlung  radiation  loss  due  to  its  higher  atomic 
number  compared  with  D.  Before  many  of  the  major  changes  to  the 
code,  the  optimum  fuel  mix  based  on  a  maximum  Isp  value  was  found 
to  be  approximately  60%  D  and  40%  3He . 

The  first  major  alteration  to  the  original  DPF  code  dealt  with 
the  usage  of  the  total  efficiency  (EFFTOT)  of  the  pinch  device. 
The  input  file  (see  Appendix  B)  for  the  program  contains  two 
efficiency  terms,  each  corresponding  to  a  fraction  of  particles 
which  are  trapped.  The  efficiency  of  the  current  sheath,  the 
snowplow,  (EFFPLW)  is  the  fraction  of  particles  which  are  swept  up 
in  the  current  sheath  discharge  as  it  runs  the  length  of  the  anode 
as  seen  in  Figure  5.  The  pinch  efficiency  (EFFPIN)  is  the  fraction 
of  particles  that  are  caught  in  the  snowplow  effect  and  then 
trapped  in  the  plasma  pinch  at  the  tip  of  the  anode.  Thus, 


EFFTOT  -  EFFPLW  x  EFFPIN  . 


(29) 


EFFTOT  is  used  in  the  calculation  of  the  particle  number  density  in 
the  pinch.  In  the  original  program,  it  was  further  included  in  the 
calculations  for  the  magnetic  induction  at  the  pinch  surface 
(BPINCH) ,  the  bremsstrahlung  power  loss,  and  the  cyclotron  power 
loss.  In  the  new  code,  the  EFFTOT  term  is  only  used  in  the 
calculation  of  the  number  of  particles  in  the  pinch. 

The  reason  for  including  EFFTOT  only  once  in  the  calculations 
is  that  its  value  propagates  through  the  program.  The  number 
density  in  the  pinch  is  used  to  calculate  BSURF  (called  BPINCH  in 
Leakeas'  program)  ,  which  in  turn  is  used  to  calculate  the  needed 
current  (IMAX)  in  the  pinch.  EFFTOT  is  also  propagated  through  the 
calculations  of  the  bremsstrahlung  and  synchrotron  radiation  power 
losses.  In  the  original  code  the  total  efficiency  was  included  in 
these  same  expressions  that  already  included  either  BPINCH  or 
NPNCH .  Since  EFFTOT  is  0.175  using  EFFPLW  =  0.70  and  EFFPIN  = 
0.25,  the  additional  inclusion  of  EFFTOT  into  an  expression  would 
decrease  the  actual  value  of  a  particular  term  by  a  factor  of  5.71. 
This  was  found  to  affect  the  performance  of  the  DPF  considerably. 

The  condition  (3  =  1,  as  mentioned  earlier  in  the  previous 
section,  formed  the  basis  for  the  next  major  change  in  the  code. 
Study  of  the  original  program  following  changes  due  to  EFFTOT' s 
essentially  multiple  inclusion  in  some  equations  showed  that  the 
calculated  BSURF  and  number  densities  resulted  in  (3  =  2.46.  This 
would  indicate  that  the  plasma  pressure  on  the  inside  of  pinch  is 
about  two  and  a  half  times  greater  than  that  of  the  magnetic 
pressure  on  the  outside.  Such  a  plasma  is  unstable,  as  it  is  not 
confined,  but  is  expanding  outward.  Thus,  the  needed  number 
densities  in  the  pinch  cannot  be  obtained. 

Therefore,  the  new  code  is  arranged  so  (3  =  1  is  a  condition 
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inherent  to  the  calculations  in  the  program.  This  is  done  in  a 
straightforward  manner.  The  initial  ion  density  in  the  pinch  is 
given  by: 


toe  P  ini  c  ^ inic 
^av^pin 


(30) 


where  nbeq  is  the  number  of  ions  in  the  pinch  in  the  beginning,  plnit 
is  the  initial  mass  density  of  the  fuel  particles  in  the  entrance 
of  the  DPF,  Vlnlt  is  the  initial  volume  of  the  focus  chamber,  mav  is 
the  average  mass  of  each  particle  within  the  pinch,  and  Vpln  is  the 
volume  of  the  pinch  itself.  From  nbeq  the  number  densities  of  the 
D  and  3He  particles  and  the  electrons  within  the  pinch  are 
determined.  Then,  using  Eq.  (2  0)  with  P  =  1,  the  magnetic 
induction  at  the  surface,  Bb,  is  found.  Bb  is  used  to  obtain  the 
value  for  IMAX  (Imax)  needed  in  order  to  create  a  magnetic  field 
strong  enough  to  contain  the  calculated  number  density  of  the  pinch 
particles.  The  result  is  that  a  much  higher  current  than 
previously  thought  necessary  is  required.  Therefore,  a  larger 
amount  of  power  is  required  of  the  capacitor  bank  in  order  to 
supply  the  needed  current.  One  finds  that  the  turbine/generator 
system  used,  with  a  turbine  efficiency  of  100%  and  a  generator 
efficiency  of  30%,  is  not  able  to  meet  this  requirement. 

Some  additional  concerns  were  also  addressed  in  writing  the 
new  code.  The  central  focus  of  the  remaining  changes  was  the 
treatment  of  synchrotron  radiation  losses.  The  original  program 
used  a  variable  CYCREFL  for  the  fraction  of  emitted  synchrotron 
energy  reflected  into  the  plasma,  assuming  a  value  of  ranging  from 
0.6  to  0.8.  As  seen  earlier  in  this  report  in  EFFECT  OF  BETA  ON 
SYNCHROTRON  RADIATION,  such  a  value  is  too  costly  for  operation  of 
the  DPF.  It  was  shown  earlier  how  much  just  letting  (3  =  0.9  can 
affect  the  losses  due  to  synchrotron  radiation  emissions.  The 
earlier  code,  with  its  (3  value  of  ~  2.46,  means  that  no  synchrotron 
radiation  could  have  been  produced.  As  seen  earlier,  no 
synchrotron  radiation  emission  occurs  for  a  |3  value  >  1.  This 
indicates  a  CYCREFL  value  of  1.00.  This  variable  is  renamed  REFSYN 
in  the  new  program.  It  is  retained  in  the  program  because  in  the 
future,  modeling  the  DPF  with  a  P  value  less  than  1.00  may  be 
desirable.  In  such  a  case,  synchrotron  radiation  would  possibly 
escape  the  plasma,  causing  one  to  be  concerned  about  how  much  of  it 
is  reflected  into  the  plasma  by  the  chamber  walls  and  nozzle 
entrance.  Until  such  a  time,  though,  the  value  of  REFSYN  has  been 
assigned  a  value  of  1.00. 

When  research  began  on  this  project,  it  was  thought  that  the 
amount  of  synchrotron  radiation  reflected  into  the  plasma  could  be 
calculated  using  the  plasma  depth  Dp  [Ref.  2] : 
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(31) 


2  _i  2 

B  2r\.2  (1  -  P)  4 
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where  rp  is  the  plasma  radius  [m] ,  7).  is  the  electrical  resistivity 
[Q-m]  of  the  reflecting  material  at  a  specified  temperature,  and 
other  variables  are  as  noted  earlier.  The  reflecting  materials 
were  the  electrodes,  assumed  to  be  made  of  copper.  After 
discovering  that  the  device  had  a  P  value  £  1,  it  was  decided  to 
forego  such  calculations,  as  p  =  1  represents  the  ideal  condition 
for  operating  this  fusion  device.  Further,  it  would  require  much 
more  time  and  detail  to  determine  how  much  of  the  synchrotron 
radiation  emissions  (which  only  occurs  for  P  <  1)  are  reflected. 
Since  the  pinch  is  fully  outside  the  DPF  device,  it  is  unlikely  any 
appreciable  emissions  would  be  reflected  to  the  pinch.  The  mixing 
chamber  (nozzle  entrance  vessel)  walls  may  be  able  to  reflect  some 
of  the  emissions,  but  this  would  only  seem  to  be  appreciable  if  the 
chamber  size  is  close  to  that  of  the  focus  device.  One  would  need 
to  study  how  reflection  is  affected  by  distance  from  the  pinch  and 
how  synchrotron  radiation  is  affected  by  transmission  through  the 
propellant  flow  leaving  the  chamber  and  entering  the  nozzle. 

One  final  area  of  concern  was  the  cooling  of  the  electrodes, 
which  are  heated  ohmically  during  the  focus  process.  It  was 
decided  that  deuterium  gas  could  be  pumped  through  a  hollow  anode 
and  possibly  the  annulus  containing  the  snowplow  between  pulses. 
The  mass  of  deuterium  necessary  to  absorb  all  of  the  Ohmic  power 
being  produced  in  a  single  pulse  is  given  by: 

-  -^=41  [kg/ a]  .  (32) 

CP.D  111 


where  Pohm  is  the  Ohmic  heating  power  [W] ,  At  is  the  length  of  time 
in  one  pulse  [s],  cp>D  is  the  specific  heat  of  deuterium,  and  AT  is 
the  change  in  temperature  [°K]  of  the  electrodes  during  one  pulse. 
Since  the  repetition  rate  used  in  the  DPF  model  is  100  times  per 
second,  this  mass  of  deuterium  is  multiplied  by  a  factor  of  100  to 
determine  mass  flow  per  second. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


This  report  has  brought  forward  the  benefit  of  having  a  dense 
plasma  focus  with  (3  =  1:  no  escape  of  emitted  synchrotron 
radiation  from  the  plasma.  In  fact,  an  appreciable  electron 
temperature  must  be  attained  before  the  characteristic  synchrotron 
radiation  reabsorption  length  exceeds  the  radius  of  the  plasma 
pinch.  It  would  be  of  possible  value  in  the  future  to  determine  if 
the  exact  shape  of  the  plasma  pinch  would  have  any  effect  on  the 
characteristic  reabsorption  length,  X0.  Further  work  in  this  area 
would  validate  whether  the  value  of  X0  in  relation  to  the  pinch 
radius  is  a  sufficient  indicator  of  whether  synchrotron  radiation 
escapes  the  plasma.  For  the  time  being,  however,  neglecting  the 
energy  losses  due  to  the  escape  of  synchrotron  radiation  is 
appropriate  for  a  P  =  1  device. 

The  synchrotron  radiation  relation  to  the  plasma  beta  is 
incorporated  into  the  dense  plasma  focus  code  listed  in  Appendix  A. 
The  new  DPF  code  is  modularized,  so  as  new  facets  are  added  to  it, 
their  incorporation  is  made  simple  through  the  insertion  of  new 
subroutines.  Slight  corrections  have  been  made  to  the  original 
code,  as  well  as  updating  it  with  new  information.  Some  test  input 
files  have  been  run,  with  the  highest  resulting  Isp  values  in  the 
range  900-1000  seconds.  The  maximum  Thrust-to-Weight  ratios  have 
been  on  the  order  of  10'4.  These  values  are  not  as  favorable  as 
those  listed  in  the  earlier  report  by  Leakeas  [Ref.  1]  .  The 
current  code  is  based  on  the  code  listing  found  at  the  back  of 
Leakeas'  report.  The  findings  in  the  earlier  report  may  have  been 
based  on  a  DPF  device  working  in  a  pulsed  manner,  rather  than  the 
near  steady-state  conditions  on  which  the  code  herein  is  based. 
Recently,  it  has  been  found  that  a  pulsed  mode  offers  higher 
performance  [Ref.  12]  .  In  the  future,  it  would  be  advisable  to 
update  the  code  to  working  in  a  pulsed  manner  and  compare  the 
results  to  the  current,  near  steady-state  device. 
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APPENDIX  A 


C  ############################################################################# 
C  #######################  A  DENSE  PLASMA  FOCUS  MODEL  ########################## 
C  ############################################################################# 

PROGRAMMER:  Larry  T.  Cox,  Jr. 

DATE:  4  JUN  1991 

LATEST  REVISION:  12  AUG  1991 


AFFILIATIONS 


School  of  Nuclear  Engineering 

Purdue  University 

West  Lafayette,  IN  47907 

Phillips  Laboratory  (AFSC)  West 
OLAC-PL/RKFE 

Edwards  AFB,  CA  93523-5000 


Based  on  the  FORTRAN  code  "THE  DENSE  PLASMA  FOCUS: 
SCENARIO"  by  Christopher  L.  Leakeas 


A  FUSION  PROPULSION 


VARIABLE  DESCRIPTIONS 


IMPLICIT  REAL* 8  (Static) 


Variable 

BOLTZ 

CPE 

CPH2 

CPHION 

FSYNAB 

GRAV 

IMAGNT 

MH2 

MUNAUT 

PI 

RHOCU 

TDISON 

TIMDIS 

WDDN 

WDDP 

WDH3 

DDN1-DDN4 

DDP1-DDP4 

DH31-DH34 


Description 

Boltzmann  constant  (1 . 381*10**-23  J/K) 

Specific  energy  of  an  electron 

"  "  of  diatomic  hydrogen 

"  "  of  a  hydrogen  atom 

Fraction  of  synchrotron  radiation  abosorbed  in  the  walls 
and  electrodes 

Acceleration  due  to  gravity  (9.80  m/s**2) 

Current  flowing  in  the  magnet  at  the  nozzle  entrance 
Atomic  mass  of  diatomic  hydrogen 
Permittivity  of  space  (4*PI*10**-7  H/m) 

Pi  (3.1415) 

Mass  density  of  copper 

Time  of  dissociation  and  ionization  of  diatomic  hydrogen 
Time  of  dissociation  for  diatomic  hydrogen 
Energy  (work)  released  in  one  D(d,n)3He  reaction 
"  "  "  in  one  D(d,p)T  reaction 

"  "  "  in  one  3He(d,p)4He  reaction 

Coefficients  for  calculation  of  DDn  rx  rate  parameter 
"  for  "  of  DDp  rx  rate  parameter 

"  for  "  of  D-3He  rx  rate  parameter 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 

c 
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c 

c 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


IMPLICIT  REAL* 8 
Variable 


BETA 

BSURF 

COEFF 

CP  CAP 

DAYS 

ECAP 

EFFTOT 

EFFPLW 

EFFPIN 

REFSYN 

REFBRE 


FLOW 

FLOWAD 

FLOWD 

FLOMOR 

FARADS 

FFBURN 

FRACD 

FRACH3 

IDOT 

IMAX 

IMAXSQ 

INDIN 

ISP 

KTPIN 

KTLOG 

LANODE 

LOOT 

LP INCH 
MASSAV 
MCAP 
MFUEL 
MFUSYS 
MMAGNT 
MPAYLD 
MPROP 
MPROST 


(Dynamic) 

Description 


Plasma  pressure  balance  coefficient  (=  1.00) 

Magnetic  induction  at  surface  of  plasma  pinch  (Tesla) 
Coefficient 

Specific  energy  of  the  capacitors 
Number  of  days  the  mission  is  to  last 
Stored  energy  in  capacitors  required  to  supply  needed 
current  IMAX 

Total  efficiency  of  focus  device  (fraction  of  injected 
particles  which  end  up  in  the  pinch) 

Efficiency  of  the  snowplow  (fraction  of  injected  particles 
caught  in  the  snowplow) 

Efficiency  of  the  pinch  (fraction  of  snowplow  particles 
trapped  in  the  pinch) 

Reflected  fraction  of  synchrotron  radiation  back  into  the 
plasma  pinch  due  to  chamber  wall  reflection 
Reflected  fraction  of  Bremsstrahlung  radiation  back  into 
the  plasma  pinch  by  means  of  seeding  the  propellant 
and/or  fuel 

Mass  flow  of  hydrogen  into  the  turbine  (kg/s) 

Mass  flow  of  hydrogen  in  excess  of  that  needed  to  generate 
the  required  power  to  charge  the  capacitor  bank  (kg/s) 

Mass  flow  of  deuterium  coolant  through  the  electrodes  (kg/s)  r 
Mass  flow  of  hydrogren  added  to  FLOW  (kg/s) 

Capacitance  of  the  capacitor  bank  (Farads) 

Fraction  of  original  D  and  3He  number  densities  consumed 
during  the  burn 

Fraction  of  deuterium  in  the  fuel  mixture  initially 
Fraction  of  Helium-3  in  the  fuel  mixture  initially 
Time  rate  of  change  of  the  capacitor  bank  current 
(Amps/second) 

Current  required  for  operation  of  the  pinch  device  (Amps) 
IMAX*  *2 

Initial  inductance  of  the  capacitor  bank  (Henries) 

Specific  impulse  of  the  DPF  system  (seconds) 

Kinetic  temperature  in  the  pinch  (keV) 

Log  base  10  of  KTPIN  used  in  rx  rate  parameter  calculations 
Length  of  the  anode  (m) 

Time  rate  of  change  of  the  capacitor  bank  inductance 
(Henries/s) 

Length  of  the  plasma  pinch  (m) 

Average  mass  of  each  atom  and/or  ion  in  the  pinch  ( 

Mass  of  the  capacitor  bank  (kg) 

Mass  of  the  D  and  3He  fuel  carried  (kg) 

Mass  of  the  pumping  and  storage  system  for  the  fuel  (kg) 

Mass  of  the  magnet  used  in  the  magnetic  nozzle  (kg) 

Mass  of  the  payload  (kg) 

Mass  of  the  hydrogen  propellant  carried  (kg) 

Mass  of  the  propellant  pumping  and  storage  system  (kg) 
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c 

c 
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c 
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MSHELD 

MTOTAL 

NBEG 

NDBEG 

NH3BEG 

NPBEG 

NTBEG 

NH4BEG 

PABSW 

PBREM 

PCAPAD 

PDELTA 


PDELTO 

PFDDN 

PFDDP 

PFDH3 

PFT 

PFTDDN 

PFTDDP 

PFTDH3 

PFTOT 

PLEFT 

PLOSS 

PLOSTO 

PMAGNT 

POHMIC 

POWELC 

POWIN 

PSYNC 

PTBREM 

PTDIS 

PTSYNC 

PWASTE 


RADPIN 


Mass  of  the  LiH  neutron  shield  (kg) 

Total  mass  of  all  DPF  system  components  (kg) 

Number  density  of  all  atoms  and/or  ions  in  the  pinch  at 
the  beginning  (#/m**3) 

Number  density  of  deuterium  atoms/ions  in  the  pinch  at  the 
beginning  (#/m**3) 

Number  density  of  Helium-3  atoms/ions  in  the  pinch  at  the 
beginning  <#/m**3) 

Number  density  of  hydrogen/protons  in  the  pinch  at  the 
beginning  (#/m**3) 

Number  density  of  tritium/tritons  in  the  pinch  at  the 
beginning  (#/m**3) 

Number  density  of  Helium-4/alphas  in  the  pinch  at  the 
beginning  (#/m**3) 

Waste  power  absorbed  in  the  walls  (MW) 

Bremsstrahlung  power  produced  (MW) 

Additional  power  needed  to  charge  capacitor  bank  (MW) 
Change  in  power  equal  to  fusion  power  produced  minus  the 
electromagnetic  radiation  power  losses  in  one  time 
interval  (MW) 

Overall  change  in  power  over  one  pulse  (MW) 

Fusion  power  produced  from  DDn  in  one  time  interval  (MW) 
Fusion  power  produced  from  DDp  in  one  time  interval  (MW) 
Fusion  power  produced  from  D-3He  in  one  time  interval  (MW) 
Total  fusion  power  produced  in  one  time  interval  (MW) 

Total  fusion  power  produced  by  DDn  in  one  pulse  (MW) 

Total  fusion  power  produced  by  DDp  in  one  pulse  (MW) 

Total  fusion  power  produced  by  D-3He  in  one  pulse  (MW) 
Total  fusion  power  produced  in  one  pulse  (MW) 

Power  remaining  after  propellant  is  fully  ionized  to  heat 
the  ionized  propellant  (MW) 

Power  lost  to  electromagnetic  radiation  during  each  time 
interval  (MW) 

Overall  power  lost  to  electromagnetic  radiation  during  the 
length  of  one  pulse  (MW) 

Power  emitted  by  the  magnet  at  the  nozzle  entrance  (MW) 
Power  absorbed  in  the  copper  electrodes  as  a  result  of 
ohmic  heating  (MW) 

Electrical  power  generated  by  the  heated  hydrogen  gas 
to  charge  the  capacitor  banks  (MW) 

Input  power  needed  by  the  focus  to  fire  (MW) 

Synchrotron  power  emitted  during  the  current  time 
interval  (MW) 

Total  bremsstrahlung  power  loss  during  one  pulse  (MW) 

Total  power  used  to  heat  and  dissociate  the  hydrogen 
propellant  (MW) 

Total  synchrotron  power  loss  during  one  pulse  (MW) 

Total  waste  power  derived  from  absorbed  bremsstrahlung, 
absorbed  sychrotron  radiation,  and  magnet  power  used 
to  heat  the  hydrogen  propellant  before  entering  the 
turbine  (MW) 

Radius  of  the  plasma  pinch  (m) 
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C23456 


RANODE 

RATREP 

RCATH 

RHOIN 

RRDDN 

RRDDP 

RRDH3 

RRDDN1 

RRDDP1 

RRDH31 

SVDDN 

SVDDP 

SVDH3 

TAV 

THICK 

THRE 

THRION 

THRNP 

THRPRO 

THRTOT 

TIMPIN 

TSTAG 

TTHROT 

VEXE 

VEXI 

VOLTS 

VTHROE 

VTHROI 

VOLINI 

VOLPIN 

XSAREA 

ZEFF 


Outer  radius  of  the  focus  anode  (m) 

Rate  of  repetition  of  the  focus  device  (pulses/s) 

Inner  radius  of  the  focus  cathode  (m) 

Initial  mass  density  of  the  fuel  at  the  entrance  to  the 
focus  device  (kg/m**3) 

Current  reaction  rate  for  D(d,n)3He  (reactions/m**3-s) 
Current  reaction  rate  for  D(d,p)T  (reactions/m**3-s) 
Current  reaction  rate  for  3He(d,p)4He  (reactions/m**3-s) 
Initial  reaction  rate  for  D(d,n)3He  (reactions/m**3-s) 
Initial  reaction  rate  for  D(d,p)T  (reactions/m**3-s) 
Initial  reaction  rate  for  3He(d,p)4He  (reactions/m**3-s) 
Reaction  rate  parameter  for  DDn  (m**3/s) 

Reaction  rate  parameter  for  DDp  (m**3/s) 

Reaction  rate  parameter  for  D-3He  (m**3/s) 

Average  temperature  of  propellant  in  the  chamber  (K) 
Thickness  of  the  LiH  neutron  shield  (m) 

Thrust  due  to  electrons  (N) 

Thrust  due  to  ions  (N) 

Thrust  due  to  non-pinch  gases  (N) 

Thrust  due  to  propellant  (N) 

Total  thrust  (N) 

Duration  of  a  pinch  (s) 

Stagnation  temperature  of  the  propellant  (K) 

Temperature  of  the  throat  (K) 

Exhaust  velocity  of  the  electrons  (m/s) 

Exhaust  velocity  of  the  ions  (m/s) 

Voltage  of  the  capacitor  bank  (V) 

Throat  velocity  of  the  electrons  (m/s) 

Throat  velocity  of  the  ions  (m/s) 

Initial  volume  of  the  fuel  in  annular  region  between  the 
anode  and  the  cathode  (m**3) 

Volume  of  the  pinch  (m**3) 

Cross-sectional  area  of  the  annular  region  between  the 
anode  and  the  cathode  (m**2) 

(Effective  atomic  number)  x  (number  of  particles  in  the 
pinch) **2 


PROGRAM  DPF 


IMPLICIT  REAL* 8  (A  -  Z) 

INTEGER  I, J, K, STEPS, RUNTYP 
INCLUDE  'BLOCK. DPF' 

DATA  NDBEG, NH3BEG, NH4BEG, NPBEG, NTBEG  /0 . 0, 0 . 0, 0 . 0, 0 . 0, 0 . 0/ 
DATA  NPIND, NPINH3, NPINH4, NPINP, NPINT  /0 .0, 0 .0, 0 .0, 0 .0, 0 .0/ 
DATA  NPINCH, NPINEL, PDELTA  /0 . 0, 0 . 0, 0 . 0/ 

DATA  PFDH3, PFDDN, PFDDP, PFT, PFTOT  /0 .0, 0 . 0, 0 . 0, 0 .0, 0 .0/ 

DATA  PFTDDN, PFTDDP, PFTDH3  / 0 . 0 , 0 . 0 , 0 . 0 / 
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DATA  PTBREM, PTSYNC, PTLOSS  / 0 . 0 , 0 . 0 , 0 . 0 / 

DATA  PDELTO, PLOSTO  /0. 0,0.0/ 

DATA  PBREM, PLOSS, PLOSS  70.0,0.0,0.0/ 

DATA  RRDDN, RRDDP, RRDH3  /0 . 0, 0 . 0, 0 . 0/ 

DATA  ZEFF  / 0 . 0 / 

CALL  INPUTDAT (STEPS, RUNTYP) 

CALL  FOCUSPAR 
CALL  REACRATE 
CALL  PULSE (STEPS) 

CALL  POWER (STEPS) 

IF  ( (POWIN+PMAGNT+ECAP) .GT.POWELC)  THEN 

PRINT  *,  'Not  enough  energy  derived  from  waste  heat  to' 
PRINT  *,  'charge  the  capacitor  bank!!!' 

ELSE 

CALL  MARSSYS (STEPS) 

CALL  RESULTS (STEPS) 

ENDIF 

STOP 

END 


BLOCK  DATA  ASSIGNMENT 


BLOCK  DATA  BLKDAT 


REAL  *  8  BOLTZ , CPE , CPH2 , CPHION, FSYNAB , GRAV, IMAGNT , MH2 , MUNAUT , P I 
REAL  *  8  RHOCU, TDISON,TIMDIS, WDDN , WDDP , WDH3 

REAL* 8  DDN1 , DDN2 , DDN3 , DDN4 , DDP 1 , DDP2 , DDP3, DDP4 , DH3 1 , DH32 , DH3 3 , DH34 


COMMON 

& 


& 

COMMON 

& 


/ CONST1 /  BOLTZ, CPE, CPH2, CPHION, FSYNAB, GRAV, IMAGNT, MH2, 
MUNAUT , P I , RHOCU , TD I SON , T IMD I S , 

WDDN, WDDP, WDH3 

/ CONST2 /  DDN1 , DDN2 , DDN3 , DDN4 , DDP 1 , DDP2 , DDP3 , DDP 4 , 

DH31 , DH32 , DH33 , DH34 


DATA  BOLTZ, CPE, CPH2, CPHION  / 1 . 3 8E-23 , 1 . 5 17E07 , 4 157 . 0 , 82 67 . 0/ 
DATA  FSYNAB, GRAV, IMAGNT, MH2  /0 . 2, 9 . 807, 3 . 18E05, 3 . 34E-27/ 

DATA  MUNAUT, PI, RHOCU  / 1 . 257E-0 6 , 3 . 1 4 15 , 1 . 673E-08/ 

DATA  TDISON, TIMDIS  /5000 . 0, 1 . 0E-07/ 

DATA  WDDN, WDDP, WDH3  /5 . 24E-13 , 6 . 4 6E-13 , 2 . 93E-12/ 

DATA  DDN1 , DDN2 , DDN3, DDN4  /0 .29811, -2 . 08296, 5.70135, -22. 08780/ 
DATA  DDP 1 , DDP 2 , DDP3 , DDP 4  /0 . 30795, -2 . 12009, 5 . 68718, -22 . 03746/ 
DATA  DH31 , DH32,DH33, DH34  /0 . 35715, -3 . 32451, 10 . 11363, -25 . 66533/ 
END 
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INPUT  DATA 


SUBROUTINE  INPUTDAT (STEPS, RUNTYP) 

IMPLICIT  REAL *8  (A  -  Z) 

INTEGER  STEPS, RUNTYP 
CHARACTER* 1 2  INFILE 
INCLUDE  ' BLOCK. DPF' 

OPEN  (UNIT=1 , FILE=' MARS . DAT' , STATUS='OLD' , FORM=' FORMATTED' , 
&  ACCESS=' SEQUENTIAL' ) 


READ 

(1, 

*) 

READ 

(1, 

*) 

REFSYN 

READ 

(1, 

*) 

READ 

(1, 

*) 

RANODE 

READ 

u, 

*) 

READ 

(1, 

*) 

RCATH 

READ 

(1, 

*) 

READ 

(1, 

*) 

LANODE 

READ 

(1, 

*) 

READ 

(1, 

*) 

RHOIN 

READ 

(1, 

*) 

READ 

(1, 

*) 

REFBRE 

READ 

(1, 

*) 

READ 

(1, 

*) 

INDIN 

READ 

(1/ 

*) 

READ 

(1, 

*) 

FARADS 

READ 

(1, 

*) 

READ 

(1, 

*) 

CPCAP 

READ 

(1, 

*) 

READ 

(1/ 

*) 

EFFPLW 

READ 

<1, 

*) 

READ 

(1, 

*) 

EFFPIN 

READ 

(1, 

*) 

READ 

(1, 

*) 

LP INCH 

READ 

(1, 

*) 

READ 

(1, 

*) 

RADPIN 

READ 

(1, 

*) 

READ 

(1, 

*) 

FRACD 

READ 

(1, 

*) 

READ 

(1, 

*) 

FRACH3 

READ 

(1, 

*) 

READ 

(1, 

*) 

KTPIN 

READ 

(1, 

*) 

READ 

(1, 

*) 

RATREP 
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READ 

(1,  *) 

READ 

(lr*) 

TIMPIN 

READ 

(lr*) 

READ 

(1,  *) 

STEPS 

READ 

(lr  *) 

READ 

(lr  *) 

DAYS 

CLOSE  (UNIT=1 ) 


RETURN 

END 


FOCUS  PARAMETER  CALCULATIONS 


SUBROUTINE  FOCUSPAR 
IMPLICIT  REAL* 8  (A  -  Z) 

INCLUDE  ' BLOCK. DPF' 

XSAREA  =  PI  *  (RCATH* *2  -  RANODE**2 )  !  [m**2] 

VOLINI  =  XSAREA  *  LANODE  ! [m**3] 

VOLPIN  =  PI  *  RADP IN* *2  *  LPINCH  !  [m**3] 

EFFTOT  =  EFFPLW  *  EFFPIN 

MASSAV  =  FRACH3*5 . 00E-27  +  FRACD*3 . 34E-27  ! [kg] 

Calculate  Initial  Pinch  Number  Density 

NBEG  =  EFFTOT  *  RHOIN  *  VOLINI  /  (MASSAV* VOLPIN)  ! [#/m**3] 

NDBEG  =  FRACD  *  NBEG  ! [#/m**3] 

NH3BEG  =  FRACH3  *  NBEG  ![#/m**3] 

Use  beta  value  of  1.0  so  that  magnetic  and  particle  pressures  balance  at  the 
plasma's  surface  and  then  calculate  the  magnetic  induction  B  at  the  surface 

BSURF  =  SQRT (  (2 . 0*NDBEG+3 . 0*NH3BEG) *1 . 6E-1 6*KTPIN*2 . 0*MUNAUT  )  ![T] 

Needed  Current  is  obtained  using  Ampere's  Law 

IMAX  =  BSURF *2 . 0  *PI*RADPIN  /  MUNAUT 
IMAXSQ  =  IMAX*  *2 

NTBEG  =0.0 
NPBEG  =0.0 
NH4BEG  =0.0 


RETURN 

END 
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REACTION  RATE  PARAMETER  CALCULATIONS 


SUBROUTINE  REACRATE 


IMPLICIT  REAL* 8  (A  -  Z) 


INCLUDE  ' BLOCK. DPF' 


KTLOG  = 

LOG10 (KTPIN 

SVDH3  = 

1.0E-06  *  1 

Sc 

DH33*KTLOG 

SVDDN  = 

1.0E-06  *  1 

Sc 

DDN3  *KTLOG 

SVDDP  = 

1.0E-06  *  1 

Sc 

RETURN 

END 

DDP3*KTLOG 

0** (DH31* (KTLOG*  *  3 ) 
+  DH34)  ! [m**3/s] 

0** (DDN1 * (KTLOG*  * 3 ) 
+  DDN4 )  ! [m**3/s] 

0** (DDP1 * (KTLOG*  * 3 ) 
+  DDP4 )  ! [m**3/s] 


+  DH32* (KTLOG*  *2 )  + 
+  DDN2* (KTLOG**2)  + 
+  DDP2* (KTLOG*  *2 )  + 


PULSE  CALCULATIONS 


SUBROUTINE  PULSE (STEPS) 

IMPLICIT  REAL  *  8  (A  -  Z) 

INTEGER  I ,  J,  K, STEPS 

INCLUDE  ' BLOCK. DPF' 

NPIND  =  NDBEG 
NPINH3  =  NH3BEG 
NPINP  =  NPBEG 
NPINH4  =  NH4BEG 

J  =  0 
PRINT  * 

DO  235  I  =  1, STEPS-1 
J  =  J  +  1 

IF  (J.EQ.500)  THEN 

PRINT  220,  'STEP  #',I 
220  FORMAT  ('+',A,I6) 

J  =  0 

END  IF 
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RRDH3 

RRDDN 

RRDDP 


NPIND  *  (NPINH3*SVDH3)  ! [m**-3*s**-l ] 

0.25  *  NPIND  *  (NPIND*SVDDN)  ! [m* *-3*s* * -1 ] 
0.25  *  NPIND  *  (NPIND*SVDDP)  ! [m**-3*S**-l ] 


Charged  particle  fusion  power  from  pinch 


& 

& 

& 


PFDH3  =  RRDH3*WDH3*VOLPIN*RATREP*TIMPIN* 
(1.0/STEPS) *1.0E-06  ! [MW] 

PFDDN  =  0 . 25*RRDDN*WDDN*VOLPIN*RATREP*TIMPIN* 
(1.0/STEPS) *1.0E-06  ! [MW] 

PFDDP  =  RRDDP *WDDP*VOLP IN* RATREP* TIMPIN* 
(1.0/STEPS) *1.0E-06  ! [MW] 

PFTOT  =  PFDH3  +  PFDDN  +  PFDDP  ! [MW] 


NPINEL  =  NPIND+2 . 0 *NP INH3+NP INT+2 . 0*NPINH4+NPINP 
ZEFF  =  (NPIND+4 . 0 *NP INH3+NP INT+4 . 0*NPINH4+NPINP) 

PSYNC  =  6.21E-17  *  VOLPIN  *  BSURF**2  *  NPINEL  * 

&  KTPIN  *  (1 . 0+KTPIN/204 . 5)  *  RATREP  *  TIMPIN  * 

&  (1.0/STEPS)  *  1.0E-06  ! [MW] 

PBREM  =  5.35E-37  *  NPINEL  *  !*  EFFTOT* *2 

&  SQRT (KTPIN)  *  VOLPIN  *  ZEFF  * 

&  RATREP  *  TIMPIN  *  (1.0/STEPS)  *  1.0E-06  ! [MW] 

PLOSS  =  (1 . 0-REFBRE) * PBREM  +  ( 1 . 0-REFSYN) *PSYNC  ! [MW] 
PDELTA  =  PFTOT  -  PLOSS  ! [MW] 

Add  the  terms  of  the  steps  to  the  power  values 


PFT  =  PFT  +  PFTOT 

!  [MW] 

PLOSTO 

= 

PLOSTO 

+ 

PLOSS 

!  [MW] 

PDELTO 

= 

PDELTO 

+ 

PDELTA 

!  [MW] 

PFTDH3 

= 

PFTDH3 

+ 

PFDH3 

!  [MW] 

PFTDDN 

= 

PFTDDN 

+ 

PFDDN 

!  [MW] 

PFTDDP 

= 

PFTDDP 

+ 

PFDDP 

S  [MW] 

PTBREM 

= 

PTBREM 

+ 

PBREM 

!  [MW] 

PTSYNC 

= 

PTSYNC 

+ 

PSYNC 

!  [MW] 

.EQ. 1) 

THEN 

RRDH31 

= 

RRDH3 

RRDDN 1 

= 

RRDDN 

RRDDP 1 

= 

RRDDP 

END  IF 


NPIND  =  NPIND  -  (RRDH3  +  2.0*RRDDN+ 

&  2 . 0*RRDDP)  *  TIMPIN  *  (1.0/STEPS)  ![#/m**3] 

NPINH3  -  NPINH3  -  (RRDH3  -  RRDDN)  *  TIMPIN  * 

&  (1.0/STEPS)  ! [#/m**3] 

NPINT  =  NPINT  +  RRDDP  *  TIMPIN  *  (1.0/STEPS)  ! [#/m**3] 

NPINH4  =  NPINH4  +  RRDH3  *  TIMPIN  *  (1.0/STEPS)  ![#/m**3] 

NPINP  =  NPINP  +  (RRDH3  +  RRDDP)  *  TIMPIN  *  (1.0/STEPS)  ![#/m**3] 

NP INCH  =  NPIND  +  NPINH3  +  NPINT  +  NPINH4  +  NPINP  ![#/m**3] 
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235  CONTINUE 


RETURN 

END 

C================================================ 

c 

C  POWER  CALCULATIONS 

C 

c================================================ 

SUBROUTINE  POWER (STEPS) 

IMPLICIT  REAL* 8  (A  -  Z) 

INTEGER  I ,  J,  K, STEPS 
INCLUDE  ' BLOCK. DPF' 


C 

C  Stored  electrical  energy  needed  in  capacitor  =  ECAP 
C 

ECAP  =  0.5  *  INDIN  *  ( IMAX+IMAGNT) * *2  *  1.0E-06  ! [MJ] 

C 

C  Determine  voltage  of  capacitor  bank  from  ECAP  value  and  FARADS  value 

C 

VOLTS  =  SQRT (2 . 0*ECAP*1 . 0E+0 6 /FARADS) 

POWIN  =  VOLTS* ( IMAGNT+IMAX) *RATREP*TIMDIS* 1 . 0E-06  ! [MW] 

C 

C  Synchrotron  radiation  power  absorbed  in  walls  and  electrodes  (in  a  BETA  = 

C  1.00  situation,  there  is  NO  synchrotron  emission,  so  PABSW  would  be  equal 
C  to  0.00  in  such  a  case) 

C 

PABSW  =  FSYNAB  *  (1.0-REFSYN)  *  PTSYNC  ! [MW] 

C 

C  Total  power  dissipated  in  electrodes  by  ohmic  heating 
C 

POHMIC  =  IMAXSQ  *  RHOCU  *  LANODE  *  (1 . 0/ (PI *RANODE**2) +1 . 0/ 

&  (PI* (RCATU**2-RANODE**2) ) )  *  1.0E-06  *  TIMDIS  *  RATREP  ! [MW] 

print  *,  ' POhmic=' , POHMIC 
C 

C  Power  removed  from  magnet 
C 

PMAGNT  =0.01  ! [MW] 

C 

C  Total  waste  power  derived  from  walls  and  magnet  (deuterium  coolant  flowing 
C  through  electrodes  escapes  out  into  the  chamber  with  the  heat  of  POHMIC 
C  in  it) 

C 

PWASTE  =  PABSW  +  PMAGNT  +  (1 . 0-REFBRE) *PTBREM  ! [MW] 

PDELTO  =  PDELTO  -  PABSW  -  POHMIC  -  PMAGNT 
C 
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C  Electric  power  is  %  of  waste  power  recovered  with  100%  efficient  turbine  and 
C  30%  efficient  generator 
C 

POWELC  =  0.30  *  PWASTE  ! [MW] 

Mass  flow  of  H2  at  2000K  into  the  turbine — minimum  flow  rate  to  avoid 
overheating  of  turbine,  if  ALL  of  PWASTE  is  needed  to  charge  the 
capacitor  and  turbine  temperature  cannot  exceed  2000K 

FLOW  =  PWASTE  *  70.0  /  (2000.0-20.0)  ! [kg/s] 

Power  needed  beyond  POWELC  for  ECAP;  1.0/0. 3  factor  is  due  to  30%  efficiency 

PCAPAD  =  ECAP  *  RATREP  -  POWELC  ! [MW] 

IF  (PCAPAD . GT . 0 . 0 )  THEN 

FLOWAD  =0.0  ! [kg/s] 

print  * , ' No  excess  propellant.' 

ELSEIF  (PCAPAD.LT. 0.0)  THEN 

FLOW  =  FLOW  +  1.0/0. 3  *  PCAPAD  *  70.0  / 

&  (2000.0-20.0)  ! [kg/s] 

FLOWAD  =  -1.0  *  1.0/0. 3  *  PCAPAD  *  70.0  / 

&  (2000.0-20.0)  ! [kg/s] 

ENDIF 

Deuterium  coolant  flows  through  anode  and  between  anode  and  cathode  to  cool 
the  electrodes  between  pulses  (assume  it  behaves  as  the  H2  gas  and  is  also 
at  2000K)  escapes  into  the  mixing  chamber. 

Amount  of  deuterium  needed  is  calculated  assuming  electrodes  are  made  of 
copper 

FLOWD  =  (1 . 6E+06*POHMIC  *  100.0  *  1.0E-04  *  0.009  *  100.0)  / 

&  (14.28  *  1000.0  *  1300.0) 

Assumed  turbine  exit  temperature  =  700K 

Gas,  having  passed  through  the  turbine  (at  700K) ,  is  injected  into  the  stream 
flowing  out  of  the  mixing  chamber,  which  includes  the  excess  flow 
of  heated  propellant  at  2000K  and  the  deuterium  coolant  flowing  out  of  the 
copper  electrodes. 

TAV  =  ( (FLOW-FLOWAD) *700.0  +  ( FLOWAD +FL0WD) *2000 . 0 )  /  FLOW  ! [K] 

Assume  gas  absorbs  all  fusion  power,  mixes  uniformly,  and  exits  at  uniform 
temperature.  Assume  gas  aborbs  heat  up  to  5000K  as  H2  with  GAMMA  =1.40  . 
Thus,  power  removed  from  system  is: 

PTDIS  =  (FLOW+FLOWD)  *  CPH2  *  (TDISON-TAV)  +  PCAPAD  ! [MW] 

IF  (PDELTO-POWIN . LE .0.0)  THEN 
PTDIS  =  0.0 

PRINT  *,  'There  is  no  power  remaining  to  heat  the  plasma' 

PRINT  *,  'gases  to  5000K.' 
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ELSEIF  (PTDIS .GE. (PDELTO-POWIN) )  THEN 
PTDIS  =  PDELTO  -  POWIN 

TAV  =  TDISON  -  (PTDIS  -  PCAPAD)  /  ( (FLOW+FLOWD) *CPH2 ) 

END  IF 

print  *,  'PTDIS (MW) =', PTDIS 

PLEFT  is  power  left  to  be  absorbed  by  an  assumed  completely  dissociated  and 
ionized  combination  of  electron  and  hydrogen  ion  gases  (at  5000K) 

IF  (TAV. GT. 5000.0)  THEN 

PLEFT  =  PDELTO  -  POWIN  -  PTDIS  ! [MW] 

ELSE 

PLEFT  =0.0 

END  IF 

PRINT  *,  'PLEFT (MW) =' , PLEFT 

Power  absorbed  by  the  gases 

IF  (TAV. LE. TDISON)  THEN 
TSTAG  =  TAV 

ELSE 

TSTAG  =  TDISON  +  PLEFT  /  ((FLOW+FLOWD)  *  (CPHION*0 . 999455  + 

&  0 . 000545*CPE) )  ! [K] 

ENDIF 

TTHROT  =  TSTAG  /  1.35  ! [K] 

VTHROI  =  SQRT (2 . 0*CPHION* (TSTAG-TTHROT) )  ![m/s] 

VTHROE  =  SQRT (2 . 0*CPE* (TSTAG-TTHROT) )  ! [m/s] 

Flow  exits  twice  as  fast  because  of  nozzle  expansion 

VEXI  =  2.0  *  VTHROI  ![m/s] 

VEXE  =  2.0  *  VTHROE  ![m/s] 

Resultant  thrust 

THRE  =  5.45E-04  *  (FLOW+FLOWD)  *  VEXE  ! [N] 

THRION  =  0.999455  *  (FLOW+FLOWD)  *  VEXI  ! [N] 

THRPRO  =  THRE  +  THRION  ! [N] 

THRTOT  =  THRPRO  ! [N] 

FFBURN  =  (NPIND+NPINH3)  /  (NDBEG+NH3BEG) 

RETURN 

END 


MARS  MISSION  &  OTHER  SYSTEM  CALCULATIONS 


SUBROUTINE  MARSSYS (STEPS) 
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IMPLICIT  REAL *8  (A  -  Z) 

INTEGER  STEPS 
INCLUDE  ' BLOCK . DPF' 

:  100  Metric  tons  payload  mass  (1  metric  ton  =  1000  kg) 

MPAYLD  =  1.0E+05  ! [kg] 

:  Mass  of  propellant  and  mass  of  propellant  system  and  structure 

MPROP  =  (FLOW+FLOWD)  *  DAYS  *  86400.0  ! [kg] 

MPROST  =  0.10  *  MPROP  ! [kg] 

:  Mass  of  deuterium  and  He-3  fuel 

MFUEL  =  8.703  *  RHOIN  *  VOLINI  *  RATREP  *  DAYS  *  86400.0  ! [kg] 

MFUSYS  =  0.1  *  MFUEL  ! [kg] 

:  Convert  ECAP  to  kJ  and  find  capacitor  mass  (noting  CPCAP  =  2  kJ/kg) 

MCAP  =  ECAP  *  1000.0  /  CPCAP  ! [kg] 

:  LiH  shield  thickness 
■* 

THICK  -  0.1  *  LOG (0.125* (RRDDN+RRDT)  *  TIMPIN  *  RATREP  *  VOLPIN 
&  *  DAYS  /  (4 .  0*PI*1 . 157E+12) )  !  [m] 

:  Density  of  LiH  =  725  kg/m**3 

MSHELD  =  725.0  *  PI  *  THICK  ! [kg] 

:  Mass  of  magnet 

MMAGNT  =67.55  ! [kg] 

:  Total  mass 

MTOTAL  =  MPAYLD  +  MCAP  +  MPROP  +  MPROST  +  MFUEL  +  MSHELD 
&  +  MFUSYS  f  MMAGNT  ! [kg] 

:  isp 

ISP  =  THRTOT  /  (GRAV  *  (RHOIN*VOLINI *RATREP+FLOW+FLOWD) )  ![s] 

RETURN 

END 

RESULTS 
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SUBROUTINE  RESULTS (STEPS) 
IMPLICIT  REAL *8  (A  -  Z) 


INTEGER  I,J,K,  STEPS 
INCLUDE  ' BLOCK. DPF' 


OPEN  (UNIT=2 , FILE=' DPF . OUT' , STATUS='NEW' ) 


WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE  (2,*) 
WRITE (2,*) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2,  *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2,  *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 


'Initial  fill  gas  density  =  ',RHOIN 

'Needed  current  (in  Amps)  =  ' , IMAX, ' ,  when' 

'  Plasma  pinch  temperature  =  ',KTPIN 

'Initial  fractions  of  D  and  3He:' 

'D  ' , FRACD/ ' 3He' , FRACH3 

'Initial  and  Final  Number  Density  (in  m**-3):' 

'D' , NDBEG, NPIND 
' 3 He ' , NH3BEG, NP INH3 
' 4He' , NH4BEG, NPINH4 
'P' , NPBEG, NPINP 
'  T'  , NTBEG, NPINT 

'Fraction  of  D-3He  orig.  #  densities  burnt  =  ',  FFBURN 

'Rx  Rate  Parameters  for  DHE3,  DDN,  DDP :  ' 

' (in  m**3/s) ' 

SVDH3 

SVDDN 

SVDDP 

'Beginning  and  ending  rx  rates  for  DHE3,  DDN,  DDP:' 

' (in  m**-3*s**-l) ' 

RRDH31 , RRDH3 
RRDDN1 , RRDDN 
RRDDP1, RRDDP 

'Charged  particle  fusion  power  of  DHE3,  DDN,  DDP:' 

'  (in  MW) ' 

PFTDH3, PFTDDN, PFTDDP 

' - >  Total  fusion  power  (in  MW)  =  ' , PFT 

'Power  needed  to  operate  focus  (in  MW):' 

POWIN 

'Power  losses  for  Brem  and  Sync  Radiation  (in  MW):' 

( 1 . 0-REFBRE) *PTBREM, (1.0-REFSYN) *PTSYNC 
'Total  radiation  power  losses  (in  MW):' 

PLOSTO 
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WRITE <2,  *) 


IF  ( 0 . 3  *PDELTO-POWIN-PMAGNT . LT .0.0)  THEN 

WRITE (2,*)  'Net  DECREASE  in  power...' 

WRITE (2,*)  PDELTO-POWIN-PMAGNT, '  [MW]' 

ELSEIF  ( 0 . 3*PDELTO-POWIN-PMAGNT . GT . 0 . 0 )  THEN 
WRITE (2,*)  'Net  INCREASE  in  power...' 

WRITE (2,*)  0 .3*PDELTO-POWIN-PMAGNT, '  [MW]' 

ELSE 

WRITE (2,*)  'NO  NET  CHANGE  in  reactor  power.' 

ENDIF 


WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2,  *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE  (2, *) 
WRITE {2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE  (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE  (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 
WRITE (2, *) 


'Total  electrical  power  produced  =  ',POWELC,'  [MW]' 

'H  mass  flow  rate  (kg/s)  =  ' , FLOW+FLOWAD 
'D  mass  flow  rate  (kg/s)  =  ' , FLOWD 

'Plasma  stagnation  temperature  (in  K,  eV) : ' 

TSTAG, TSTAG/12000 . 0 

'Ion  and  electron  exit  velocities  (in  m/s):' 

VEXI, VEXE 

'Final  propellant  thrust  (in  N) : ' 


THRPRO 

'Firing  time  =  ',DAYS,'  days' 

' ======================The  System=====================' 

'  Area  Mass  (kg)  ' 

/  _  _ _ _  / 

' PAYLOAD  ' , MPAYLD 

' PROPELLANT  ' , MPROP 

'PROPELLANT  SYSTEM  AND  STRUCTURE' , MPROST 
'  FUEL  ' , MFUEL 

'FUEL  SYSTEM  ',MFUSYS 

'CAPACITOR  ' , MCAP 

'SHIELD  ' , MSHELD 

' MAGNET  ' , MMAGNT 

f _  t 

'TOTAL  MASS  ' , MTOTAL 

t - f 


'Total  thrust  (in  Newtons,  Lbf):' 

THRTOT, THRTOT/4 .4482 

'Thrust  to  Weight  Ratio  =', THRTOT  /  (MTOTAL *GRAV) 

' Isp  (Specific  Impulse)  =  ',ISP,'  [ s ] ' 

'Thrust  *  Isp  =  ', THRTOT  *  ISP,'  [N-s]' 

'!!!!!!!! END !!!!!!! OF !!!!!!! ! OUTPUT !!!!!!! DATA !!!!!!!!' 


CLOSE  (UNIT=2 ) 
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RETURN 

END 


APPENDIX  B 


'NCHROTRON  RADIATION  REFLECTION  COEFFICIENT 

00 

IODE  RADIUS  IN  METERS  (0.0508) 

0508 

iTHODE  RADIUS  IN  METERS  (0.0800) 

1050 

IODE  LENGTH  IN  METERS  (0.382) 

382 

IITIAL  FILL  GAS  DENSITY  IN  KG/M~3  (2.2E-4) 

2E-4 

IFLECTED  PBR  FRACTION  DUE  TO  SEEDING  OF  PLASMA  FUEL 
00 

IITIAL  INDUCTANCE 
2E-8 

UPACITOR  BANK  SPECIFIC  ENERGY  IN  KJ/KG 
00 

JOWPLOW  EFFICIENCY  FACTOR 
70 

:nch  efficiency  factor 
25 

:nch  length  in  meters 
0254 

:nch  radius  in  meters 
5E-3 

(ACTION  OF  DEUTERIUM  IN  INITIAL  FILL  GAS 
55 

(ACTION  OF  HELIUM-3  IN  INITIAL  FILL  GAS 
45 

(NTED  OPERATING  TEMPERATURE  OF  PINCH  (keV) 

).0 

jASMA  focus  firing  rate  in  S**-l 

)0.0 

JRATION  OF  STABLE  PINCH  IN  SECONDS 
.  00E-4 

JMBER  OF  ITERATIONS  (>=  1000  for  accurate  modeling) 

)00 

)TAL  FIRING  TIME  (for  entire  mission  to  Mars  w/delta-v  of  30  km/s)  IN  DAYS 

)0 .0 
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